Characterizing quark gluon plasma by thermal photons and lepton pairs 
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The photon spectra measured by the ALICE collaboration in Pb+Pb collisions at Large Hadron 
Collider (LHC) energies has been analyzed with a view of extracting the properties of thermal 
system formed in these collisions. The results of the analysis are compared with the previously 
studied spectra measured at Super Proton Synchrotron (SPS) and Relativistic Heavy Ion Collider 
(RHIC) energies. The thermal dilepton spectra from the Pb+Pb collision at LHC energy has been 
predicted for the initial conditions constrained by the thermal photon spectra at the same collision 
conditions. The slope of the photon has been used to estimate the rise in the effective degeneracy 
with the charged particle multiplicity of the system. The slopes of the lepton pair spectra for 
different invariant mass windows have been used to conjecture the average radial flow velocity of 
the high temperature phase. 

PACS numbers: 25.75.-q,25.75.Dw,24.85.+p 



I. INTRODUCTION 

The main aim of colliding heavy nuclei at relativis- 
tic energies is to create a hot and/or dense thermal 
phase of matter where the quarks and gluons are not 
confined inside hadrons but move within nuclear vol- 
ume. Such a phase of matter is called quark gluon 
plasma (QGP). The detection of QGP in heavy ion 
collisions (HIC) at RHIC and LHC energies is one of 
the most challenging jobs both for experimentalists 
and theorists working in this field primarily because 
of the extremely transient nature of the QGP. The 
QGP evolves dynamically in space and time due to 
high internal pressure, due to which the system cools 
and reverts to hadronic matter. The electromag- 
netically interacting particles, photons and dileptons 
PHI] (see for review) are considered to be pene- 
trating probes of the matter formed in the collisions 
of heavy ions at relativistic energies because (i) they 
are produced at each space-time point of the system 
and (ii) their mean-free paths are much larger than 
the system-size and hence bring the information of 
the production points very efficiently. 

The efficiency of photons for being considered as 
an competent probe of QGP largely depends on the 
ability to disentangle the photons produced from 
various stages of evolution of the system formed in 
HIC. Therefore, first we identify the possible sources 
of photons above those coming from the decays of 7r° 
and rj mesons etc, as provided by the data because 
photons from these decays are already eliminated 
and hence need not be considered in the present 
analysis. Photons produced from the evolving mat- 
ter under consideration are: 

(i) due to the initial hard collisions of the partons 
from the nucleons of the colliding nuclei. This con- 
tribution may be estimated by using the techniques 
of perturbative QCD (pQCD) and the data from pp 
collisions may be used to validate such calculations. 



The px distributions of photons from proton+proton 
(pp) collisions at a given energy can be used as a 
benchmark for the hard contribution in HIC. There- 
fore, estimation of these contributions with minimal 
model dependence is important. In view of this, in 
the present analysis we estimate the hight px con- 
tributions in HIC by using the experimental data 
obtained in pp collisions |7] as described in section 
III. 

(ii) Photons are also produced from the interac- 
tions of the 'yet-to-be equilibrated' partons i.e. from 
the time span between the collision point and the on- 
set of thermalization. In case where the thermaliza- 
tion time scale is very small the contributions from 
this interval will be insignificant and hence can be 
neglected. 

(iii) Thermal photons originating from the interac- 
tions of the (a) quarks and gluons in the bath and (b) 
thermal hadrons (7r, p, rj, u>, a\ etc). The estimation 
of the thermal contribution depends on the space- 
time evolution scenario that one considers. In case 
of a deconfinement phase transition,which seems to 
be plausible at RHIC energies (see Q for a review) , 
one assumes that QGP is formed initially. The equi- 
librated plasma then expands, cools, and reverts to 
hadronic matter and finally freezes out at a tempera- 
ture, Tf ~, called freeze-out temperature. Evidently 
there will be thermal radiation from QGP as well as 
from the luminous hadronic fireball which has to be 
estimated as accurately as possible in order to have 
a reliable estimate of the initial temperature. 

The momentum distributions of photons (and 
dileptons) produced from a thermal system depend 
on the temperature (T) of the source through the 
thermal phase space factors of the participants of the 
reactions [9(- Consequently, the transverse momen- 
tum (jpr) spectra of photon reflects the temperature 
of the source. For an expanding system the situation 
is, however, far more complex. The thermal phase 
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space factor changes - by several factors e.g. the 
transverse kick received by low px photons due to 
flow originating from the low temperature hadronic 
phase (realized when T < T c ) populates the high px 
part of the spectra • As a consequence the inter- 
mediate or the high px part of the spectra contains 
contributions from both QGP as well as hadrons. 
The photon spectra measured experimentally repre- 
sents the space-time integrated yield from the matter 
that evolves from an initial hot and dense phase to a 
comparatively cooler and diluted phase of hadronic 
gas. Therefore, the temperature extracted from such 
spectra will exhibit the average temperature of the 
system. 

The experimental data on the px distributions of 
photons from various collision energies and colliding 
systems have been analyzed by using different kinds 
of models [IiT - [l8j mainly to extract the temperature 
of the hot phase formed after the collisions. In the 
present work the slope of the photon spectra has 
been used to estimate the enhancements of the ef- 
fective degeneracy with increased multiplicity. The 
slope of dilepton spectra for different invariant mass 
windows have been used to extract the average ra- 
dial flow velocity of the high temperature phase. 

In the next section we very briefly describe the 
various mechanisms of photon and lepton pair pro- 
ductions. Section II is dedicated a compressed de- 
scription of the space-time dynamics of the system. 
The results are presented in section III and section 
IV is devoted to summary and discussions. 



0(aa s ) has been done by resuming ladder diagrams 
in the effective theory [211 ] . which has been used in 
the present work. A set of hadronic reactions with all 
possible isospin combinations have been considered 
for the production of photons [22T-I24} from hadronic 
matter. The effect of hadronic dipole form factors 
has been taken into account in the present work as 




B. Sources of thermal dileptons 

The dominant source of the thermal dileptons 
from QGP is the qq annihilation (25|. For the low 
mass dilepton production from HM the decays of 
thermal light vector mesons namely p, u and <f> have 
been considered. The change of spectral function of 
p due to its interaction with n, ui, a%, hi (see [26|, [27| 
for details) and baryons [28[ have been included in 
evaluation of lepton pairs from HM. For the spec- 
tral function u) the width at non-zero temperature 
is taken from Ref. [29( and no medium effect has 
been considered for <f>. The continuum part of the 
spectral function of p and lj have also been included 
in the dilepton production rate [f| H3] • The model 
employed in the present work leads to a good agree- 
ment with NA60 dilepton data [3l[ for SPS collision 
conditions [32| . 

II. EXPANSION DYNAMICS 



A. Sources of thermal photons 

The emission of thermal photons and lepton pairs 
from QGP and thermal hadrons have been discussed 
in earlier works. In the present work we briefly out- 
line various processes for the emissions and refer to 
the literature for the details. The rate of thermal 
photon production per unit space-time volume per 
unit four momentum volume is given by (see 
H for a review): 
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d 3 p (2tt) 3 



ImIL*J BE (E,T) 



(1) 



where ImlT^ is the imaginary part of the retarded 
photon self energy and fsE(E,T) is the thermal 
phase space distribution for Bosons. For an expand- 
ing system, the energy E should replaced by u^p^, 
where p M and are the four momentum and the 
hydrodynamic four velocity respectively. 

The Hard Thermal Loop [l!| approximations has 
been used by several authors [20] to evaluate the 
photon spectra originating from the interactions of 
thermal quarks and gluons. The complete calcula- 
tion of emission rate of photons from QGP to order 



The space time evolution of the system formed in 
Pb+Pb collisions at ^Jsnn — 2.76 TeV has been 
studied by using relativistic hydrodynamics with 
longitudinal boost invariance [33J and cylindrical 
symmetry [34| . We assume that the system reaches 
the state of equilibrium at a time Ti after the colli- 
sion. The initial temperature, Ti can be related to 
the measured hadronic multiplicity (dN /dy) by the 
following relation for system undergoing isentropic 
expansion: 



^ = 7rR\4a q TS/c 



(2) 



where Ra is the radius of colliding nuclei, c is a 
constant ~ 4 and a q = (tt 2 /90)g q where g q (= 2 x 
8 + 7x2x2x3x Np/8) is the degeneracy of quarks 
and gluons in QGP, -/Vp^number of flavours. The 
value of dN/dy can be calculated from the following 
equation [351 ]: 



dN dn P p < N part 

dy ~ [ X) dy 2 
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< N coU > 

Ncoii is the number of collisions and contribute to 
x fraction to the multiplicity dn pp /dy measured in 
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TABLE I: The values of various parameters - thermal- 
ization time (ri), initial temperature (Tj) and hadronic 
multiplicity dN/dy - used in the present calculations. 

ysjviv 2.76 TeV 

centrality 0-40% 
^ 1212 

dy 

n 0.1 fm 

Ti 553 MeV 

T c 175 MeV 

Tf 100 MeV 

EoS Lattice QCD 



TABLE II: The values the (pr) for different collision 
energies obtained from the direct photon data at low 
Pt by fitting with ao x exp[— pr/ai]. 



(l/N part )(dN ch /dr,) 
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FIG. 1: Transverse momentum spectra of direct photon 
at 2.76 TeV energy for Pb+Pb collision at 0-40% cen- 
trality. 



pp collision. The number of participants, N par t con- 



tributes a fraction (1 



of dn p p/dy. The values 



of N par t and N co u are estimated by using Glauber 
Model and the results are in agreement with (36j . 
We have used dn pp /dy = 4.31 and x = 0.1 at 
y?SNN =2.76 TeV. It should be mentioned here that 
the values of dN/dy (through N part and N co u in 
Eq. 13) and hence the T (through dN/dy in Eq. [2) 
depend on the centrality of the collisions. The val- 
ues of Ra for a given centrality has been evaluated 

1/3 

by using the relation - Ra ~ l.l N' r + . 

We use the lattice QCD EoS |37[ for the QGP 
phase and hadronic resonance gas EoS for the 
hadronic phase [38j . The kinetic freeze out tempera- 
ture, Tf = 100 MeV is constrained by the px spectra 
of hadrons. The ratios of various hadrons measured 
experimentally at different ^/snn indicate that the 
system formed in heavy ion collisions chemically de- 
couple at T c h(> Tf). Therefore, the system remains 
out of chemical equilibrium from T cri to Tf. The 
deviation of the system from the chemical equilib- 
rium is taken in to account by introducing chemical 
potential for each hadronic species [39(. 



III. RESULTS AND DISCUSSION 



A. pr distributions of photons and dileptons 
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FIG. 2: The variation of (pr) with the increase in mul- 
tiplicity for different collision energies. 
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The direct photon spectra from Pb+Pb collisions 
is measured at ^/sjvjv = 2.76 TeV. However, no 
data at this collision energy is available for pp in- 
teractions. Therefore, prompt photons from p+p 



FIG. 3: The invariant mass distribution of thermal lep- 
ton pairs for Pb+Pb collision at LHC. 
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FIG. 4: The transverse mass distribution of thermal lep- 
ton pairs for Pb+Pb collision at LHC. 
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FIG. 5: The variation of effective slope with invariant 
mass bins for Pb+Pb collision at LHC. The dashed line 
is obtained by setting radial velocity, vt = 



collision at ^/snn — 7 TeV has been used to esti- 
mate the hard contributions for nuclear collisions at 
^/snjv = 2.76 TeV by using the scaling (with ^/snn) 
procedure used in |7j. For the Pb+Pb collisions the 
result has been scaled up by the number of collisions 
at this energy (this is shown in Fig.[T]as prompt pho- 
tons) . The high px part of the data is reproduced by 
the prompt contributions reasonably well. At low px 
the hard contributions underestimate the data indi- 
cating the presence of a possible thermal source. 

The thermal photons with initial temperature ~ 
553 MeV along with the prompt contributions ex- 
plain the data well (Fig. [IJ, with the inclusion of 
non-zero chemical potentials for all hadronic species 
considered [39[ (see also (ioj). 

It is well known that transverse momentum spec- 
tra of photons act as a thermometer of the interior of 
the plasma. The inverse slope of the thermal distri- 
bution is a measure of the average (over evolution) 
effective (containing flow) temperature of the sys- 
tem. We have extracted the average effective tem- 



perature (~< pt >) from the thermal distributions 
of photons at different collision energies - i.e. for 
SPS, RHIC and LHC energies. Fig[2]shows the vari- 
ation of (pt) with multiplicity for different collision 
energies. To minimize the centrality dependence of 
the results the dN c h/drj is normalized by N part . The 
results clearly indicate a significant rise in the aver- 
age pt (< Pt >) while going from SPS to RHIC to 
LHC. The values of (pt) for different collision en- 
ergies are given in the table [TTJ Since photons are 
emitted from each space time point of the system, 
therefore, the measured slope of the pt spectra rep- 
resents the average effective temperature of the sys- 
tem. 

The quantity, P°//(= l/N part dN ch /dr)) is pro- 
portional to the entropy density. Therefore, 
p"jy/< pt > 3 oc <?eff, the average effective statisti- 
cal degeneracy, a quantity which changes drastically 
if the colour degrees of freedoms deconfmed i.e. if a 
phase transition takes place in the system. We find 
that the entropy density (s ~ g e ffT 3 ) at LHC in- 
creases by almost 96% compared to RHIC and there 
is an enhancement of 46% at RHIC compared to 
SPS. However, part of this increase is due to the in- 
crease in the temperature and part is due to increase 
in degeneracy. To estimate the increase in the degen- 
eracy we normalize the quantity p^Jf by < pt > 3 - 
Therefore, we estimate < Pt > 3 from the anal- 
ysis of the experimental data and found that there 
is a 15% increase in this quantity from SPS to RHIC 
and 35% increase from RHIC to LHC. 

We evaluate the invariant mass (M) spectra, and 
the transverse mass (Mt) spectra of lepton pairs 
with initial conditions, EoS etc constrained by the 
measured photon spectra at ^/JWn — 2.76 TeV. 
The emission processes from QGP and the hadronic 
phases are taken from Ref. |26[, therefore, we do 
not repeat the details here to save space. The 
Mt(= \/ M% v + p't) spectra is evaluated for differ- 
ent M windows (M ranging from Mi to M%, with 
M av = (Mi + M 2 )/2). The M spectra displayed in 
Fig. |3] indicates that by selecting Mi and Mi ap- 
propriately, one can extract the properties of QGP 
(M av > 1.5 GeV) or hadronic system (M av ~ m p , p 
mass). Therefore, for example, the slope of the Mt 
spectra (Fig. gj at M av ~ 2 GeV and 0.77 GeV pro- 
vide information about the average temperature and 
flow of the QGP and hadronic phases respectively. 

In Fig. [5] the slope of the Mt spectra, T e t f for 
different M windows has been plotted. To under- 
stand the effect of flow on the slope of the spectra 
we evaluate the spectra and hence estimate the slope 
by switching on and off the radial flow. The results 
with (solid line) and without (dashed line) flow are 
depicted in Fig. [5] The slope, T e /y with flow (v r ^ 0) 
may be parametrized as T e tf = Tf + M av < v r > 2 . 
The difference in the slope due to non-zero v r and 
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the observation of the dominance of the QGP phase 
at large M help in estimating the radial flow of the 
QGP phase. The estimated value of < v r >~ 0.065 
for M = 2 GeV. A very small value of 0.065 is jus- 
tified because the M ~ 2 GeV range correspond to 
very early time when flow is not fully developed. 
However, for M ~ 1.6 GeV we found < v r >~ 0.24 
where QGP phase contributions dominate, indicat- 
ing the fact that the QGP formed at LHC collision 
conditions undergo significant radial flow. Similarly, 
the value of < v r > for the hadronic phase (near the 
p peak) is 0.46. The value of < v r > in the hadronic 
phase at freeze-out will be more than 0.46 as this 
is the average value of the hadronic matter. The 
< iv >~ 0.37 for M ~ 1.3 GeV where the QGP and 
the hadronic matter contribute almost equally (see 
Fig.©. 

IV. SUMMARY AND DISCUSSIONS 

In summary, we have analyzed the photon spectra 
measured by ALICE collaboration in Pb+Pb colli- 
sions at y^JJvjv = 2.76 TeV with a view of extracting 



the properties of thermal system formed in these col- 
lisions. The slope of the photon spectra have been 
compared with the previously measured spectra at 
Super Proton Synchrotron (SPS) and Relativistic 
Heavy Ion Collider (RHIC) energies. The thermal 
lcpton pair spectra from the Pb+Pb collision at LHC 
energy has been estimated with the same initial con- 
dition that is used to reproduce the thermal photon 
spectra at the same collision conditions. The in- 
crease in the effective statistical degeneracy at RHIC 
and LHC relative to SPS have been estimated from 
the slope of the photon spectra. The radial flow ve- 
locity of the QGP and hadronic phases have been 
assessed from the invariant mass and transverse mo- 
mentum distributions of the lepton pairs. As the 
data from ALICE collaboration is well reproduced 
by the sources of photons described above, photons 
from other sources e. g. due to jet-thermal parton 
interactions [4l| and induced emissions by the hard 
partons due to multiple interactions in the QGP [42[ 
are ignored in the present analysis. These processes 
do not appear to be essential for the present range 
of transverse momentum. 
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